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We develop an ensemble density functional theory for the fractional quantum Hall effect using a lo- 
cal density approximation. Model calculations for edge reconstructions of a spin-polarized quantum 
dot give results in good agreement with semiclassical and Hartree-Fock calculations, and with small 
system numerical diagonalizations. This establishes the usefulness of density functional theory to 
study the fractional quantum Hall effect, which opens up the possibility of studying inhomegeneous 
systems with many more electrons than has heretofore been possible. (Contact mdj@physics.ucf.edu. 



Cond-mat paper cond-mat/9506141 .) 
The fractional quantum Hall effect (FQHE) is manifested in a two-dimensional electron gas (2DEG) in a strong 
ly-) • magnetic field perpendicular to the plane of the electrons [Q. The effect is due to the electron-electron interactions, 
Q\ \ which cause a downward cusp in ground state energy as a function of filling factor v — 2it£ 2 B n at certain rational 
fillings v = p/q. Here £b — \JUc/(eB) is the magnetic length, B the magnetic field strength, and n the electron 
density. These downward cusps give rise to an energy gap, and so for these values of v, the ground state of the system 
Oh 1 is an incompressible liquid. Experimentally, FQHE systems are accessible in strip geometries, Corobino geometries, 
and more recently in quantum dots with as few as ~ 50 electrons, as well as wide- well hetero junctions and 
double layer systems. Theoretically, different aspects of FQHE systems can be modeled by Laughlin's wavefunction 
by Hartree-Fock or composite Fermion Hartree by semiclassical methods ||^-[l0[|, by field theoretical 

approaches El]], which also exist for the edge excitations on the boundary of FQHE systems OL and by exact 
numerical diagonalizations. At the present, numerical diagonalizations are limited to systems with of the order of 10 
electrons. It is highly desirable to have a computational approach which accurately treats inhomogeneous systems 
with the order of 10 2 -10 3 electrons. One such approach which is in principle valid for any interacting electron system 
is the density functional theory (DFT) [|l3|-^5|. 

We have developed an ensemble DFT scheme within the local density approximation (LDA) for the fractional 
quantum Hall effect, and applied it to spin-polarized circularly symmetric quantum dots [jf?]]. The results are in good 
agreement with results obtained by semiclassical [p|~ p0|, Hartree-Fock |JJ|] (for cases where the correlations do not 
play a major role), or exact diagonalization methods |16|]. Our calculations show that the exchange and correlation 
effects are very well represented by the LDA and that our approach provides a computational scheme to model large 
inhomogeneous FQHE systems. We note that there exist previous formal DFTs for strongly correlated systems, in 
particular for high-temperature superconductors jl8| , and DFT calculations of high-T c materials |l!| and transition- 
metal oxides [^0|. Fcrconi, Gellcr and Vignalc |l(| have also recently studied FQHE systems within the DFT using an 
extended Thomas-Fermi approximation, including a LDA for the exchange-correlation energy. However, ours are, to 
the best of our knowledge, the first practical DFT-LDA calculations of a strongly correlated system in strong magnetic 
fields, and demonstrate the usefulness of the DFT-LDA in studying large inhomogeneous FQHE systems. 

In typical DFT calculations of systems of -/V e i electrons, the standard Kohn-Sham (KS) scheme (2lJ is implemented, 
in which the particle density n(r) is expressed in terms of a Slater determinant of N > N e \ KS orbitals, ip a (jr). 
These obey an effective single-particle Schrodinger equation H c gip a = e a ip a , which is solved self-consistently by 
occupying the N c \ KS orbitals with the lowest eigenvalues e a (we identify the Fermi energy of the system with 
the largest e a of the occupied orbitals), and iterating. This scheme works well in practice for systems which are 
noninteracting u-representable jl4{l5|, i.e., systems for which the true electron density can be represented by a single 



Slater determinant of single-particle wavefunctions. However, when the KS orbitals are degenerate at the Fermi 
energy there is an ambiguity in how to occupy these degenerate orbitals. This is the case for the FQHE, as we 
now demonstrate. Consider an FQHE system in the xy-pl&ne with the magnetic field along the z-axis. A circularly 
symmetric external potential V r cx t(r) = VextC?") (due, e.g., to a uniform positive background charge density) confines 
the systems such that the density is fixed with a local filling of v = 1/3 up to an edge at r$ (tq ^> £b) where the 
density falls to zero within a distance of order £b- That such systems exist is well demonstrated by the excellent 
agreement between the Laughlin wavefunction and experiments, and by many numerical calculations p2| , ^6[ . Due to 
the circular symmetry we can label single-particle orbitals by angular momentum to, and by a "band" or Landau level 
index n > 0. The orbitals ^ m , n {v) are centered on circles of radii r m w \j2m£B with Gaussian fall-offs for r <C r m 
and r 3> r m . The single-particle orbitals with n = are then in the bulk all degenerate, and the degeneracy is not 
lifted by electron-electron interactions since the system is homogeneous in the bulk. In order to obtain a constant 
density at v = 1/3 even at the center of the system, all single- particle orbitals in the bulk with n = must have 
occupancies 1/3. If the Fermi energy lay above the energies of the bulk orbitals, they would all be filled and one 
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would have v = 1. Therefore, to get occupancies 1/3 the Fermi energy must lie at the degenerate energy e m o of these 
orbitals. Thus, in applying DFT to the FQHE we can expect a huge degeneracy of KS orbitals at the Fermi energy. 
Consequently, the particle density cannot be expressed in terms of a single Slater determinant. Instead, the density 
has to be constructed from an ensemble of Slater determinants, i.e., the orbitals at the Fermi energy are assigned 
fractional occupation numbers, just as the Laughlin wavefunction is not a single Slater determinant, but a highly 
correlated state with average occupancies of 1/3 of single-particle states. This is generally known as ensemble density 
functional theory ]l^ , ^5[ . 

Although ensemble DFT has been developed formally, there are in practice few examples of applications and 
calculations using ensemble DFT for ground state calculations. A significant aspect of our work is that we have 
developed an ensemble scheme which is practical and useful for the study of the FQHE. In ensemble DFT, any 
physical density n(r) can be represented by n(r) = J2mn /mn|'0mn(r)| 2 , where f mn are occupation numbers satisfying 
< fmn < 1, and the orbitals ip mn satisfy the equation 
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A(r) + Vext(r) + Vk(r) + ^ xc (r,B) V ip {m . n) {v) = e m „Vw(r), (1) 



where V x A(r) = B(r). In equation (Q), Vn(r) is the Hartree interaction of the 2D electrons, and V KC (r, B) is the 
exchange-correlation potential, defined as a functional derivative of the exchange-correlation energy E xc [n(r), B] of the 

system with respect to density: V^ c (r, B) = ^^ffeffi 2 ^ ■ We will hereafter not explicitly indicate the parametric 
dependence of V xc and E xc on B. The question is then how to determine these orbitals and their occupancies in 
the presence of degeneracies. Here, we have devised a scheme to obtain a set of occupancies which (a) converges 
to physical densities (to the best of our knowledge) for FQHE systems, and (b) reproduces finite temperature DFT 
as well as the standard KS scheme for noninteracting w-rcpresentable systems. This scheme may be of much more 
general applicability to general systems which are not noninteracting f-representable other than the FQHE. In our 
scheme, we start with input occupancies and single-particle orbitals and iterate the system N eq times using the KS 
scheme. The number N eq is chosen large enough (about 20-30 in practical calculations) that the density is close to 
the final density after N eq iterations. Were the system noninteracting i>-representable, we would now essentially be 
done. However, in this system there are now in general many degenerate or near-degenerate orbitals at the Fermi 
energy, and small fluctuations in the density between iterations cause the KS scheme to occupy a different subset of 
these orbitals each iteration. This corresponds to constructing different Slater determinants each iteration. While 
the occupation numbers f mn of these orbital are zero or unity more or less at random each iteration, the average 
occupancies, i.e., the occupancies averaged over many iterations, become well defined and approach the value, say, 1/3 
for orbitals localized in a region where the local filling factor is close to v = 1/3. We use this to construct an ensemble 
by accumulating running average occupancies f mn after the initial N eq iterations and use these to calculate densities. 
Thus, our algorithm essentially picks a different (near-)degenerate Slater determinant after each iteration, and these 
determinants are all weighted equally in the ensemble. It is clear that this scheme reproduces the results of the KS 
scheme for noninteracting w-representable systems (for which the KS scheme picks only the one Slater determinant 
which gives the ground state density) for 7V eq large enough. We have numerically verified that a finite-temperature 
version of our scheme converges to a thermal ensemble at finite temperatures down to temperatures of the order 
of 10~ 3 TiOj c /kB- We have also performed preliminary Monte Carlo simulations about the ensemble obtained by our 
scheme. The results are that to within numerical accuracy our scheme gives the lowest energy. 

In the LDA, the exchange-correlation energy is assumed to be a local function of density, E KC /N = J dre xc (i/)n(r), 
where e xc (v) is the exchange-correlation energy per particle in a homogeneous system of constant density n = v /(2ir£ B ) 
and filling factor v. Experience has shown that the LDA often works surprisingly well, even for systems in which the 
electron density is strongly inhomogeneous (l|. In the FQHE, the length scale of exchange-correlation interactions 
and density fluctuations is given by the magnetic length £b due to the Gaussian fall-off of any single-particle basis in 
which the interacting Hamiltonian is expanded. The densities are relatively smooth on this length scale, which gives 
us additional hope that the LDA will work well for the FQHE, too. For the exchange-correlation energy per particle 
of a uniform electron gas in a constant magnetic field we use the Pade approximant |^3j e xc (^) = exc ^ + " + e *£ ^ ntyU ^ | 
where ej c c is the zero- magnetic field result [ p4| . The term e xc (^) consists of two terms. The first one is a smooth 
interpolation formula pH] e™ M {v) between ground state energies at some rational fillings. The second one, e xc (i / ), 
is all-important for the study of the FQHE. This term contains the cusps in the ground state energy which cause the 
FQHE. Here we have used a simple model which captures the essential physics. We model e^ c {v) by constructing it to 
be zero at values of v = p/q which display the FQHE. Near v = p/q, e xc (^) is linear and has at v = p/q a discontinuity 
in the slope related to the chemical potential gap A/i = q(\A p \ + |A^|). Here A p ,h are the quasiparticle (hole) creation 



energies which can be obtained from the literature 26 27j at fractions v = p/q. Farther away from v = p/q, e xc (^ 
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decays to zero. Finally, in the LDA V xc (r) is obtained from t xc {v) as V xc (r) — d["<gc(")] ^ cons t an t B. In our 
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calculations, we restrict ourselves to include only the cusps at v — 1/3, 2/5, 3/5 and v — 2/3, which are the strongest 
fractions. 

A technical difficulty arises in the LDA: The discontinuities in V xc (r) in the LDA give rise to a numerical instability. 
The reason is that an arbitrarily small fluctuation in charge density close to an FQHE fraction gives rise to a finite 
change in energy. To overcome this problem, we made the compressibility of the system finite, but very small, 
corresponding to a finite, but very large, curvature instead of a point-like cusp in e xc at the FQHE fractions. What 
we found worked very well in practice was to have the discontinuity in chemical potential occur over an interval of 
filling factor of magnitude 10~ 3 . This corresponds to a sound velocity of about 10 6 m/s in the electron gas, which 
is three orders of magnitude larger than the Fermi velocity of a 2D electron gas at densities typical for the FQHE. 
Figure [l] depicts a plot of V xc as a function of filling factor used in our calculations. 

We have self-consistently solved the KS equations Eqs. (|l|) for a spin-polarized quantum dot in a parabolic external 
potential, V cx t(r) = ^m*fl 2 r 2 by expanding the KS orbitals Vw(r) = e™ 1 ^ ip mn (r) in the eigenstates of H = 

(p — |A(r)) 2 in the cylindrical gauge, A(r) = ^Br4>, including the four lowest Landau levels (n = 0, . . . , 3). We 
chose eo — 13.6, appropriate for GaAs, and a confining potential of strength Q htt = 1.6 meV. In particular, we 
have used our DFT-LDA scheme to study the edge reconstruction of the quantum dot as a function of magnetic field 
strength. As is known from Hartree-Fock and exact diagonalizations P,[q,|l6|,|6"|,|r| , for strong confinement the quantum 
dot forms a maximum density droplet in which the density is uniform at v = 1 in the interior, and falls off rapidly to 
zero at r w \J2~NIb = Tq. As the magnetic field strength increases, a "lump" of density breaks off, leaving a "hole" 
or deficit at about r = ro. This effect is due to the short-ranged attractive exchange interaction: it is energetically 
favorable to have a lump of density break off so that the system can take advantage of the exchange energy in the 
lump. As B is further increased, the correlations will cause incompressible strips with densities v = p/q to appear 
[p|,p| p^JTo[ | on the edges, and incompressible droplets to form in the bulk at densities v = p/q. Figure || depicts various 
stages of edge reconstruction obtained by us as the magnetic field strength is increased. The value of fl for which the 
exchange lump appears compares very well with the value found by De Chamon and Wen [g] in Hartree-Fock and 
numerical diagonalizations. At higher fields still, the incompressible strips appear at the edges, and incompressible 
droplets are formed in the bulk. We emphasize that incompressible regions that appear in our calculations are not 
due to the presence of a uniform positive background density which tends to fix the bulk density at the value of the 
background density. 

There is also another edge effect caused by correlations. For particular, stiff confining potentials, so-called composite 
edges J2^,[l(| can appear. These can be thought of as particle-hole conjugates of uniform incompressible droplets. 
Consider a droplet with a bulk density corresponding to v = 1/3, falling off to zero at the edge. A incompressible 
droplet with a bulk density of v = 2/3 is obtained by particle- hole conjugation. However, at the edge, the density will 
first rise to v — 1 (since the density of the v — 1/3 droplet drops to zero), and then eventually drop to zero. Note 
that this argument is based on particle-hole conjugation, which is an exact symmetry of the lowest Landau level [p0[ , 
and it is unclear if composite edges exist in real systems, which do not strictly obey particle-hole symmetry. 

Figure || depicts the occupations of the Kohn-Sham orbitals and the particle density (inset) for a system where the 
confining potential is supplied by a uniform positive background charge density n + = 2/(6 , k£%) for r < ro, where ro 
is fixed by charge neutrality. From this figure, we see that for this choice of potential, the system forms a composite 
edge, even though our system does not obey particle-hole symmetry. We therefore conclude that such edges can exist 
in real systems. We have also verified the stability of all our incompressible regions by adding or subtracting a particle 
from the system. 
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FIG. 1. Exchange-correlation potential V xc as function of filling factor in units of e 2 /(eo^s) for < v < 1. The increase in 
V^c at an FQHE filling occurs over a range of filling factor of 0.004. 

FIG. 2. Edge reconstruction of a quantum dot as the magnetic field strength is increased. Plotted here is the local filling 
factor v(r) for a parabolic quantum dot with HQ =1.6 meV and 40 electrons. For magnetic field strengths B < 2.5 T the dot 
forms a maximum density droplet, and for B ~ 3.0 T, an exchange hole is formed. For stronger magnetic fields, incompressible 
regions form, separated by compressible strips. 

FIG. 3. Occupancies f m o us orbital centers r m for a composite edge of a system of 45 electrons (diamonds). Here B = 5.0 
T. Near the edge, the occupancies rise to unity. The full line shows the corresponding local filling factor. 
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